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1.  Introduction  and  Problem  Statement 

In  work  described  below,  we  have  investigated  the  scattering  properties  of  one¬ 
dimensional,  randomly  rough  surfaces  that  produce  multiple  scattering.  It  has  been  shown 
that  a  surface  with  reasonable  statistical  properties  (such  as  Gaussian  statistics  and  a  Gaussian 
correlation  function)  may  produce  backscattering  enhancement  and  related  polarization  effects. 
The  purpose  of  this  project  has  been  to  study  the  scattering  properties  of  such  surfaces  using 
rigorous  theoretical  methods  which  take  into  account  multiple  scattering  and  the  finite 
conductivity  of  the  surface.  In  the  calculations  to  be  discussed  below,  surface  parameters  are 
chosen  so  as  to  provide  a  direct  comparison  with  experimental  results,  thus  providing  a  direct 
means  of  determining  the  current  state  of  agreement  or  disagreement  between  theory  and 
experiment. 


2.  Theoretical  Discussion 

To  second  order  in  the  scattered  amplitudes,  the  properties  of  any  scatterer  may  be 
completely  specified  by  the  matrix  relating  the  incident  and  the  scattered  Stokes  vectors,  which 
we  will  refer  to  as  the  Stokes  matrix  of  the  scatterer.  For  a  planar  surface,  rough  in  one 
dimension,  that  is  illuminated  in  the  plane  perpendicular  to  the  grooves  of  the  surface,  it  is 
well-known  ti.at  the  Stokes  matrix  5  takes  on  the  simplified  form 

* 11  *12  0  0 

5  m  *12  *n  0  0 

0  0  J33  su 

.  0  0  -J34  J33  . 


so  there  are  only  four  unique  elements  that  are  of  interest  These  distinct  matrix  elements  are 
defined  as 


b 


s„=  [<l/lll2>  +  <l/22l2>)/2  , 

*, 2=  [<!flll2>~<!/22l2>]/2  . 

^33  =  [  </ll  fli  ^  1  » 

and  j34  =  -  Im  [  </n  f22  >  ]  , 

where /n  and  fn  are,  respectively,  die  p-  and  s-polarized  scattered  amplitudes,  and  the  angle 
brackets  denote  an  a\  erage  over  the  ensemble  of  surfaces.  All  matrix  elements  are  implicitly 
functions  of  the  incident  angle  60  and  of  the  scattering  angle  6S. 

The  Stokes  matrix  elements  have  been  calculated  for  rough  gold  surfaces  at  two  different 
wavelengths  (X  *  1.152  and  3.392  p.m).  The  optical  properties  of  gold  at  these  wavelengths 
were  described  by  the  complex  refractive  indexes  (n+ik)  =0.312+i  7.93  and  1.96+i  20.7 , 
respectively.  The  statistical  properties  of  the  surface  profiles  were  similar  to  those  measured  in 
two  recently  fabricated  surfaces.  The  stochastic  profile  had  Gaussian  statistics  and  a  standard 
deviation  of  the  surface  heights  denoted  a;  the  height  correlation  function  was  Gaussian  with  a 
transverse  correlation  length  a.. 

The  approach  used  is  based  on  the  numerical  calculation  of  expressions  for  the  scattered 
field  obtained  from  the  exact  field  equations  and  from  approximate  boundary  conditions.  This 
Impedance  Boundary  Condition  method  has  been  shown  to  give  reasonable  results  for  the 
surface  statistics  and  refractive  indices  considered.  The  incident  beam  was  of  Gaussian  form 
with  transverse  width  2g  .  A  realization  of  the  surface  profile  having  the  prescribed  statistical 
properties  was  first  evaluated  at  Nx  (=500)  equally  spaced  points,  the  values  of  the  incident 
and  scattered  field  at  the  surface  were  then  computed  at  these  points,  and  the  scattered  field  was 
obtained  as  a  function  of  the  scattering  angle  for  this  realization  of  the  surface  profile.  This 
procedure  was  repeated  for  Np  (=4000)  realizations  of  the  surface  profile,  and  the  realizations 
of  the  scattered  field  obtained  in  this  fashion  were  used  in  the  determination  of  the  second 
moments  of  the  scattered  field  and  of  the  Stokes  matrix  elements.  The  number  of  points  Nx 


and  the  total  width  of  the  surface  used  in  the  calculations  were  increased  until  a  sufficient 
accuracy  (1%)  was  obtained  on  all  four  matrix  elements;  this  typically  occurred  for  klAx  -10 
(where  Ac  *  L!  Nx)  and  L/g=4.  The  ratio  of  the  total  scattered  power  divided  by  the  total 
incident  power  was  0.98  and  0.96  for  a  s  and  a  p-polarized  incident  beam,  respectively. 

3.  Results 

In  Fig.  1  we  show  the  matrix  elements  Jn(0s),  *i2(e»>.  r33(0s),  and  r34(0s)  calculated  in 
the  cases  0o=O°,  10°,  and  30°,  for  the  surface  71  IN,  and  for  an  incident  wavelength  of  1.152 
pm.  Gold's  refractive  index  is  (n  +  i  k)  *  (.312  +  7.93  0  at  this  frequency.  The  ratio  of  the  He 
correlation  radius  a  to  the  incident  wavelength  X  is  3.1,  and  the  ratio  of  the  standard  deviation 
of  the  surface  heights  o  to  a  is  0.55.  The  surface  roughness  is  sufficiently  strong  to 
extinguish  completely  the  specular  part  of  the  scattered  field.  The  matrix  element  ,  displays  a 
peak  in  the  retroreflection  direction  0S«  0o,  and  this  peak  as  a  decreasing  magnitude  with 
respect  to  the  background  of  diffuse  intensities  when  the  angle  of  incidence  increases.  For 
small  angles  of  incidence  this  backscattering  peak  has  two  secondary  maxima  located  on  either 
side  of  the  peak,  at  about  15°  from  the  retroreflection  direction.  The  presence  of  a 
backscattering  peak  and  the  other  features  observed  in  this  matrix  element  have  been  previously 
reported  in  calculations  of  the  s-  and  p-polarized  cross-sections  of  similar  surfaces.  Since  the 
matrix  element  su  is  the  sum  of  these  cross-sections,  the  observed  behavior  may  be  explained 

in  terms  of  single  and  multiple  scattering  in  a  similar  fashion  as  in  the  case  of  these  polarized 
cross-sections.  The  element  r,2  is  about  an  order  of  magnitude  smaller  than  r,,,  but  it  is 

significant  and  it's  angular  dependence  emerges  clearly  out  of  the  statistical  noise  still  present 
after  averaging  over  Np  realizations  of  the  surface  profile.  This  element  is  the  difference 
between  the  previously  considered  r-  and  p-polarized  cross-sections,  and  its  angular 
dependence  includes  seve;al  sign  changes.  The  matrix  element  r33  has  a  magnitude  which  is 
similar  to  the  magnitude  of  the  other  diagonal  elements  off  the  Stokes  matrix  (r, ,),  and  it  also 
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Figure  1.  Theoretical  calculations  of  the  Stokes  matrix  elements  for  a  rough  surface  of 
roughness  o  *  1.95  pm,  1/e  correlation  radius  3.57  pm,  illumination  wavelength  X  *  1.152 
pm,  and  gold’s  refractive  index  (n  +  t  <)  =  (.312  +  7.93  i) .  They  are  plotted  as  a  function  of 
scattering  angle  for  normal  incidence  (top),  10°  incidence  (center),  and  30°  incidence  (baiom). 
Circles  denote  Ju,  squares  denote  r» 2,  diamonds  denote  *33,  and  triangles  denote  J34. 
Compare  with  Figure  2. 
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Figure  2.  (For  Comparison  with  Figure  1)  Experimental  results  for  the  Stokes  matrix  elements 
for  gold  surface  71  IN.  Parameters:  Roughness  o  *  1.95  pm.  Me  correlation  radius  3.57  pm, 
illumination  wavelength  X  ■  1.152  pm.  They  are  plotted  as  a  function  of  scattering  angle  for 
normal  incidence  (top),  10°  incidence  (renter),  and  30°  incidence  (bottom).  Circles  denote  Sn. 
squares  denote  sjj,  diamonds  denote  J33,  and  triangles  denote  S34.  Normalization  assumes 
unit  power  of  the  scattered  s-polarized  intensity. 


displays  a  very  similar  structure  in  the  backscattering  direction.  This  matrix  element  is  negative 
almost  everywhere  except  in  the  retroreflecdon  direction  for  small  angle  of  incidence.  This 
behavior  may  be  understood  if  one  consider  that  gold  is  highly  conducting  at  this  frequency.  It 
was  demonstrated  in  our  earlier  work  that,  for  a  perfectly  conducting  surface  with  similar 
statistical  parameters,  the  main  contributions  to  s33  are  single  and  double  scattering 
contributions  (with  shadowing),  as  is  the  case  for  rn,  but  that  in  this  matrix  elements  the 
single  scattering  contributions  add  negatively  while  the  double  scattering  contributions  add 
positively.  The  results  shown  in  Fig.  1,  demonstrate  thus  the  separability  of  the  single  and 
double  scattering  contributions  in  the  case  of  highly  conducting  surfaces.  The  effects  of  the 
finite  conductivity  of  the  surface  may  be  clearly  observed  in  the  matrix  element  r34.  Even 

though  the  angular  dependance  of  this  matrix  element  are  similar  when  computed  for  perfectly 
and  finitely  conducting  surfaces,  the  magnitude  of  su  as  obtained  for  the  gold  surface  at  this 

frequency  is  more  than  twice  that  obtained  in  the  perfectly  c  nducting  case. 

In  Fig.  2  the  experimental  data  obtained  for  the  surface  71  IN,  and  for  an  i.icident 
wavelength  of  1.152  pm,  are  presented  in  the  cases  0o*O°,  10°.  and  30°;  these  data  may 

thus  be  compared  to  the  theoretical  results  presented  in  Fig.  1.  It  is  seen  that  the  main 
discrepancies  observed  in  this  cases  are  the  following:  the  results  for  r,,  indicate  that  the 

measured  scattered  intensity  is  systematically  larger  at  high  scattering  angle  than  the  calculated 
one,  and  related  differences  occur  in  the  other  matrix  elements;  the  enhanced  backscattering 
peak  is  slightly  narrower  in  the  experimental  data  than  in  the  theoretical  results,  and  the  strength 
of  the  multiple  scattering  contributions  relative  to  the  single  scattering  contributions  is  larger  in 
the  experimental  data.  Some  of  these  discrepancies  will  be  addressed  in  Figs.  5  and  6. 

The  theoretical  results  obtained  in  the  case  of  the  the  surface  7 1 1 N  illuminated  at  a 
wavelength  X  «  3.392  pm  are  presented  in  Fig.  3.  As  in  the  case  presented  in  ^ig.  1  a 

backscattering  peak  located  about  the  retrorefkction  direction  is  present  in  the  matrix  elements 
s,,  and  Sjj  .  The  magnitude  of  these  peaks  is  seen  to  decrease  when  the  angle  of  incidence  is 

increased.  The  width  of  the  enhancement  around  0*«  0ois  significantly  larger  for  the  longer 


0.8 


Figure  3.  Theoretical  calculations  of  the  Stokes  matrix  elements  for  a  rough  surface  of 
roughness  a  ■  1.9S  |im,  1/e  correlation  radius  3-57  pm,  illumination  wavelength  X  »  3.392 
pm.  and  gold’s  refractive  index  (it  ♦  l  k)  » (1.958  ♦  20.7  0 .  They  are  plowed  as  a  function  of 
scattering  angle  for  normal  incidence  (top),  ’.0°  incidence  (center),  and  30°  incidence  (bottom). 
Circles  denote  *u.  squares  denote  *12.  diamonds  denote  J33.  and  triangles  denote  J34. 
Compare  with  Figure  4. 
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Figure  4.  (For  Comparison  with  Figure  3)  Experimental  results  for  the  Stokes  matrix  elements 
for  gold  surface  71 1 N.  Parameters:  Roughness  a  =  1.95  pm,  1/e  correlation  radius  3.57  pm, 
illumination  wavelength  X  =  3.392  pm.  They  are  plotted  as  a  function  of  scattering  angle  for 
normal  incidence  (top),  10°  incidence  (center),  and  30°  incidence  (bottom).  Circles  denote  s\\, 
squares  denote  J12,  diamonds  denote  S33,  and  triangles  denote  S34.  Normalization  assumes 
unit  power  of  the  scattered  r-polarized  intensity. 


wavelength;  this  is  consistent  with  the  interpretation  of  the  enhancement  as  occurring  from  the 
double  reflection  of  the  light  within  the  grooves  of  the  surface:  the  typical  diffraction  width  of 
such  a  scattering  process  is  X/a .  The  off-diagonal  elements  are  z.  ~  *  an  order  of  magnitude 

smaller  than  the  diagonal  elements  but  their  angular  dependence  emerges  clearly  from  the 
statistical  fluctuations.  Contrary  to  the  case  presented  in  Fig.  1,  s12  has  a  narrow 

backscattering  peak  ( the  vertical  scale  used  in  Fig.  3  does  not  allow  for  the  resolution  of  this 
peak).  The  magnitude  and  the  angular  deper.  ience  of  the  element  s34  are  different  when 

computed  for  the  finitely  conducting  surface  and  for  the  perfectly  conducting  surface.  Fig.  4 
shows  the  experimental  data  corresponding  to  the  theoretical  results  of  Fig.  3  for  the  sake  of 
comparison.  The  same  kind  of  discrepancies  observed  when  comparing  Figs.  1  and  2  are 
occurring  in  this  case.  Th  mean  scattered  intensities  measured  at  high  scattering  angle  are 
systematically  larger  than  the  calculated  intensities,  and  the  multiple  scattering  contributions  to 
the  mean  intensities  are  stronger  in  the  experimental  data  than  in  the  calculations. 

Next  we  present  calculations  of  the  matrix  elements  performed  for  surface  parameters 
corresponding  to  surface  226J.  This  surface  is  characterized  by  a=3.5  pm  and  o=1.7  pm,  and 
is  thus  statistically  quite  similar  to  surface  71  IN.  These  calculations  were  performed  because  a 
slightly  different  method  was  used  in  its  fabrication,  resulting  in  a  higher  quality  surface  with 
an  actual  profile  closer  to  the  one-dimensional  assumption  used  in  the  calculations.  It  may  be 
seen  by  comparing  the  theoretical  results  for  this  surface  (Fig.  5)  to  the  experimental  ones 
(Fig.  6)  that  some  of  the  discrepancies  have  been  reduced.  In  particular  the  mean  intensities 
measured  at  high  scattering  angles  for  X  =  1.152  pm  are  now  very  close  to  their  calculated 
counterpan.  There  are  still  important  differences  regarding  the  strength  of  the  multiple 
scattering  contributions  to  the  mean  intensities. 

The  results  presented  in  Fig.  7  describe  the  Stokes  matrix  occurring  in  the  scattering  of  a 
beam  of  light  with  a  wavelength  X  »  3.392  pm.  These  theoretical  results  resembles  the  ones 
obtained  for  the  surface  71  IN  at  this  wavelength  (Fig.  3),  except  that  there  are  slightly  smaller 


Figure  S.  Theoretical  calculations  of  the  Stokes  matrix  elements  for  a  rough  surface  of 
roughness  0  *  1.7  urn,  1/e  correlation  radius  3.5  nm,  illumination  wavelength  X  *  1.152  Jim, 
and  gold's  refractive  index  (n  +  i  k)  «  (.312  +  7.93  i) .  They  are  plotted  as  a  function  of 
scattering  angle  for  normal  incidence  (top),  10°  incidence  (center),  and  30°  incidence  (bottom). 
Circles  denote  s\\,  squares  denote  *12.  diamonds  denote  J33,  and  triangles  denote  534. 
Compare  with  Figure  6. 
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Figure  6.  (For  Comparison  with  Figure  5)  Experimental  results  for  the  Stokes  matrix  elements 
for  gold  surface  226J.  Parameters:  Roughness  a  -  1.7  pm,  1/e  correlation  radius  3.5  Jim, 
illumination  wavelength  X  =  1.152  pm.  They  are  plotted  as  a  function  of  scattering  angle  for 
normal  incidence  (top),  10°  incidence  (center),  and  30°  incidence  (bottom).  Circles  denote  jj  j, 
squares  denote  $12,  diamonds  denote  J33,  and  triangles  denote  J34.  Normalization  is  from  first 
principles. 


Figure  7.  Theoretical  calculations  of  the  Stokes  matrix  elements  for  a  rough  surface  of 
roughness  o  ■  1.7  jun,  1/e  correlation  radius  3.5  Jim,  illumination  wavelength  X  *  3.392  Jim. 
and  gold's  refractive  index  (n  +  i  k)  ■  (1.958  +  20.7  i)  .  They  are  plotted  as  a  function  of 
scattering  angle  for  normal  incidence  (top),  10®  incidence  (center),  and  30°  incidence  (bottom). 
Circles  denote  £u>  squares  denote  J12,  diamonds  denote  *33,  and  triangles  denote  *34. 


multiple  scattering  contributions  due  to  a  smaller  ratio  a/a.  At  present  time  there  are  no 
experimental  data  available  for  comparison,  but  experiment  will  be  performed  to  that  effect. 


4.  Conclusions 

In  this  work  we  demonstrated  that  theoretical  results  may  be  obtained  for  the  Stokes  matrix 
of  a  one-dimensional  randomly  rough  surface  producing  enhanced  backscattering  when  this 
surface  is  finitely  conducting.  The  results  of  these  calculations  are  found  to  be  in  good 
agreement  with  experimental  data,  even  though  some  discrepancies  regarding  the  amount  of 
multiple  scattering  occurring  at  the  surface  have  been  determined.  Moreover  these  results 
confirm  the  significance  of  the  Stokes  matrix  elements  of  such  surfaces.  The  comparisons  of 
these  results  with  the  corresponding  calculations  obtained  in  the  case  of  perfectly  conducting 
surfaces  allow  us  to  evaluate  the  role  of  the  finite  conductivity  of  the  surface,  and  to  study 
multiple  scattering  effects  in  the  scattering  of  light  from  highly  conducting  surfaces.  For 
illumination  wavelengths  in  the  infrared,  the  impedance  boundary  condition  was  used  for  the 
first  time  in  the  kind  of  simulation  calculation  that  are  performed  in  the  study  of  enhanced 
backscattering. 
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